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ABSTRACT 

We present analysis of MACHO Alert 95-30, a dramatic gravitational microlensing 
event towards the Galactic bulge whose peak magnification departs significantly from 
the standard point-source microlensing model. Alert 95-30 was observed in real-time 
by the Global Microlensing Alert Network (GMAN), which obtained densely sampled 
photometric and spectroscopic data throughout the event. We interpret the light-curve 
"fine structure" as indicating transit of the lens across the extended face of the source 
star. This signifies resolution of a star several kpc distant. 

We find a lens angular impact parameter 6 m i n /9 source = 0.715 ± 0.003. This 
information, along with the radius and distance of the source, provides an additional 
constraint on the lensing system. Spectroscopic and photometric data indicate 
the source is an M4 III star of radius 61 ± 12Rq, located on the far side of the 
bulge at ~ 9 kpc. We derive a lens angular velocity, relative to the source, of 
21.5 ± 4.9kms _1 kpc -1 , where the error is dominated by uncertainty in the source 
radius. Likelihood analysis yields a median lens mass of O.67^g' 46 M0, located with 
80% probability in the Galactic bulge at a distance of 6.93^ 2 ' 2 5 kpc. If the lens is a 
main-sequence star, we can include constraints on the lens luminosity. This modifies 
our estimates to Mi cns = O.53^ O ' 35 M0 and Aens = 6.57l 2 ' 25 kp c - 

Spectra taken during the event show that the absorption line equivalent widths of 
Ha and the TiO bands near 6700A vary, as predicted for microlensing of an extended 
source. This is most likely due to center-to-limb variation in the stellar spectral lines. 
The observed spectral changes further support our microlensing interpretation. These 
data demonstrate the feasibility of using microlensing limb crossings as a tool to probe 
stellar atmospheres directly. 

Subject headings: dark matter - gravitational lensing - stars: low-mass, brown dwarfs - 
stars: late-type - stars: atmospheres 



1. Introduction 

The MACHO collaboration is undertaking an extensive search for gravitational microlensing 
by objects in the Galactic halo, bulge, and disk. Nightly observations of millions of stars in the 



Large Magellanic Cloud (LMC) and Galactic bulge have yielded a total of 8 LMC ( Alcock et al 



1996c]) and more than 120 bulge events ([Alcock et al. 1995a| ; |Alcock et al. 1996b| ; MACHO 



Alert system [|). Similar observations are being undertaken towards the LMC by EROS (Renault 



et al. 1996] ) with 2 reported events, as well as towards the bulge by OGLE ( Paczyhski et al. 



1995) with 18 reported events and DUO ( [Alard et al. 199~5| ) with some 10 reported events. New 



efforts to detect the microlensing of unresolved stars in M31 are being undertaken by the AGAPE 
( Gondolo et al. 1996J ) and VATT/Columbia QCrotts 1996| ) collaborations. Reviews of the field of 
gravitational microlensing are presented in |Gould (1996| ) and Paczyhski (1996). 



Statistical analysis of an ensemble of microlensing events provides a useful discriminant 
between Galactic structure models. Two years of observations towards the LMC yield a 
microlensing optical depth representing approximately 50% of a "standard" Galactic halo 



comprised entirely of massive compact halo objects ( Alcock et al. 1996c ). Similar analysis of 
Galactic bulge data, as reported by MACHO and OGLE ( [Alcock et al. 1996b ; Udalski et al. 
1994a| ), indicate an optical depth a factor of ~ 3 larger than predicted by axisymmetric Galactic 



models (e.g. priest et al. 199~l| ; Paczyhski 1991). Consistent Galactic models have recently been 
constructed which include a Galactic bar viewed nearly end-on ( Kiraga &: Paczyhski 1994 ; [Zhac| 
et al. 1996; Binney et al. 1996] ). However, it is difficult to identify clearly a lensing population 
because the mass, velocity, and distance of each lens are not uniquely determined in the standard 
microlensing solution. 

Gravitational microlensing is characterized by the transient, achromatic brightening of a 
background star due to gravitational deflection of its light by a massive "lens" passing near 
our line of sight to the source. This results in the distortion of the source disk into multiple 
(unresolved) images whose total brightness is greater than that of the original source brightness. 
Given the assumptions of a point-mass lens, point source, and unaccelerated motion, the source 
magnification A{t) is of the simple form 



A(t) 



u 2 + 2 



a 



\Ju 2 + 4 ' 



(1) 



where u(t) = ^ u 2 min + ( 2(i - t )/t ) 2 . 



These equations describe the event's "light -curve" as a function of time. A(t) is the observed 
source magnification, t = 2Re/v± is the characteristic event timescale, where v± is the transverse 



1 Current information on the MACHO Collaboration's Alert events is maintained at the WWW site 



http://darkstar.astro.washington.edu. 



velocity of the lens, and 



R E 
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(2) 



is a conveniently defined distance scale in the lens plane (the lens's Einstein radius). M is the 
mass of the deflector, and x is the ratio of the observer-lens to observer-source (L) distances. The 
equation for u(t) describes the impact parameter of the lens, scaled by Re, as it passes near the 
source, and to an d u m [ n are the values of t and u at peak magnification. 

Fitting a normal microlensing event provides only one parameter, t, containing information 
about the lens. However, t is a function of the three unknowns M, x and v±. (Strictly speaking, 
the source distance L is also unknown, but the uncertainty in L is typically much smaller than 
that in x, so it is simpler to treat L as known). This indicates a continuum of event parameters 
which can conspire to produce similar duration microlensing events. 



2. Light Curve Fine Structure and Extended Source Effects 

In instances where the standard approximations break down, light-curve "fine structure" 
describes deviations from the standard model whose nature may break this degeneracy. Important 



microlensing events in this regime include the detection of binary microlensing (Udalski et al. 



1994b|; Pratt et al. 1995|; |Alard, Mao, & Guibert 1995|; |Alcock et al. 1996d|), and the observation 



of parallax in a gravitational microlensing event ( Alcock et al. 1995b| ). An exciting possibility is 
the detection of short timescale deviations in an event light -curve due to the presence of planets 
in the "lensing zone" of the microlensing system ( pVIao fc Paczyhski 1991 ; Gould fc Loeb 1992 ; 



Bennett & Rhie 199€ 



The effects of extended source size may also become apparent in the limit of a large source, 
high magnification event, where the lens impact parameter is of the order of the projected source 



radius ( |Gould 1994| ; |Nemiroff k Wickramasinghe 1994] ; |Witt fc Mao 1994| ; [Peng 19971 ) . In this 
situation, an extra parameter is included in the microlensing solution, u* = xR*/Re- Here i?* is 
the radius of the source star, and u* describes its projection into the lens plane, scaled by Re- 
Due to its geometric nature, u(t) is not changed by this parameterization. An analytic solution 



for A{t) in this limit is presented in Witt &: Mao (1994). 



Fitting the light -curve with the extra parameter u* provides a second constraint equation for 
the three lens parameters - the angular size of the source in terms of the lens's Einstein radius. If 
we estimate the linear radius and distance of the source star from photometry and spectroscopy, 
we have an estimate of the lens proper motion relative to the source, which leads to a unique 
mass-distance relation for the lens. Such a situation is rare, and would be expected in only ~ 5% 
of events towards the galactic bulge ( Gould fc Welch 1996| ). 



To account for limb-darkening of the source, an appropriate limb-darkening law and 
coefficients must be determined for the source star. The magnification is then integrated over 
the face of the star, properly weighted by the brightness profile. The limb-darkening law we will 
choose is of the form (e.g. Claret, Diaz-Cordoves, &: Gimenez (1995] )) 



/ A (/,)// A (l) = 1 - o(l - m) - &(1 " M) 2 - (3) 

Here fi is the cosine of the angle between the observer's line of sight and the emerging stellar 
radiation, and a and b are model parameters dependent primarily upon the effective temperature 
and surface gravity of the star. 



3. Global Microlensing Alert Network (GMAN) 

To provide both the temporal and photometric resolution necessary to distinguish light -curve 
"fine structure", coordinated follow-up observations are being undertaken through the Global 
Microlensing Alert Network (GMAN) ( [Pratt et al. 1995fc |Alcock et al. 1996a| ), as well as the 



PLANET collaboration ( |Albrow et al. 19971 ). 



GMAN utilizes telescope resources at Cerro Tololo Inter- American Observatory Q (CTIO), 
Mount John University Observatory (MJUO), University of Toronto Southern Observatory 
(UTSO), and Wise Observatory which respond nightly to the dynamic status of ongoing 
microlensing events (see Table [l]). Observations from MJUO are carried out by the MOA 
collaboration of Japan and New Zealand, which is undertaking a general survey of microlensing 
events towards the Magellanic Clouds and the Galactic bulge ( |Abe et al. 1997 ). The observational 



flexibility available at these locations provided nearly continuous coverage of this event. 

Microlensing follow-up images are automatically processed on site. Photometry is carried 
out using IRAF (] scripts which call Daophotll ( Stetson 1994a ). Reductions are usually finished 



within 6 hours of image acquisition. Team members have immediate network access to GMAN 
data, which were subsequently used for real-time scheduling during this event. Normalization is 
performed on-the-fly using a list of reference stars obtained from the MACHO data set. The 
reported Daophotll error bars for UTSO and WISE photometry have been increased by a factor 
of 2.0, to account for the typical scatter in a time series of constant stars from these images. The 
MJUO data used in this analysis were obtained manually from the raw frames after the event, 
using Daophotll operating within the IRAF environment; the uncertainties were estimated by 



2 Operated by the Association of Universities for Research in Astronomy, NOAO, under cooperative agreement 
with the National Science Foundation 

3 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under contract to the National Science Foundation. 



examining the variations of standard star ratios within and between frames. Data from CTIO 
were re-reduced using Airframe (Stetson 1994b). 



The system and data collection of the MACHO experiment are described in Alcock et al. 
(1996c). The data are obtained simultaneously in two specially designed wide-pass filters with 
effective wavelengths of approximately 5150 A and 6900 A. Transformations from the MACHO 
instrumental photometry system to the standard V and Rkc system have been derived by 
comparisons with several thousand tertiary standard stars calibrated in the LMC. Checks of our 
photometry with published photometry in Baade's Window confirm a zero point uncertainty of 
less than 0.1 mag in V and R, where this uncertainty is a conservative estimate of small zero point 
differences between fields monitored by MACHO in the bulge. 



4. Photometry of MACHO Alert 95 30 

The source star in Alert 95-30 is located at a =18:07:04.26, 5 = -27:22:06 (J2000). The 
MACHO project's identifier for this star is 101.21821.128. Alert 95-30 was detected by the 
MACHO Alert system at A ~ 1.8 on Jul 24, 1995, approximately 22 days before the observed 
peak. 

The location of the source star in an optical color-magnitude diagram (CMD), V = 16.21 and 



(V — R) = 1.39, indicates this star passes the "clump giant" cuts defined in Alcock et al. (1996b ). 
Figure [l] shows a single epoch CMD for a 2x2 arcmin field surrounding the source star, which is 
indicated with a circle. We may be reasonably sure then that this star is located in the Galactic 
bulge and has a large radius. Microlensing fits to the rising portion of the light-curve predicted 
a high-magnification event for the giant source, which presented the possibility of source star 
resolution. 

Subsequent observations indicated deviations from the standard microlensing model at 
A ~ 20, approximately 2 days before the projected peak. At this time, we mounted an aggressive 
program of photometry and spectroscopy to study these deviations thoroughly. Follow-up 
observations continued nightly past the observed peak until the star returned to its baseline state. 
Intermittent observations followed to determine accurately the baseline flux as measured at each 
site. 

Figure ^ displays the MACHO 1995 bulge season light-curve of Alert 95-30 []. The alert 
date is indicated with an arrow, after which time the light-curve becomes heavily sampled with 
follow-up data. Data from all follow-up observatories have been included after determining 
the baseline flux in each passband. No microlensing fits are included in the figure, but it is 



4 We report an uncatalogued CCD trap which allowed several contaminated MACHO-red observations to pass 
through our processing stream. The majority of these were removed from the data set based upon their proximity to 
the bad columns and psf F WHM. 



apparent that the data conform to the symmetric and achromatic shape expected of gravitational 
microlensing events, at least for magnification A < 20. 

This conclusion is confirmed after fitting a normal microlensing curve to the combined 
dataset. Event parameters for this fit (fit 1 hereinafter) are listed in Table |2j, and \ 2 statistics for 
each passband are listed in Table 0. However, near the peak of the event the data deviate from the 
expected light-curve (see Fig 0). This type of deviation is not unexpected in a high magnification 
microlensing event. We note that the standard point-source microlensing model allows infinite 
magnification when the source and lens are aligned. When the extended size of the source is 
considered, the magnification is limited, as the entirety of the source disk cannot be perfectly 
aligned with the lens. If we include an extra parameter in the fit to account for the angular size 
of a constant surface brightness source star (fit 2), we reduce \ 2 by 1084. The deviations near 
the peak of MACHO 95-30 are significantly reduced with this model, which we interpret as clear 
justification for the extended source microlensing interpretation. 

The data in Table [3| indicate an extended source fit \ 2 P er degree of freedom of approximately 
1 for the GMAN follow-up data. The implication here is that the data scatter around the fit is 
what one would statistically expect due to measurement error. Thus, the data are in excellent 
agreement with the extended source microlensing fit, considering our photometry is at the 1 — 2% 
level, as indicated in column 3 of Table |3|. However, the MACHO data exhibit x\of between 
3 — 4, which formally indicates a poor fit to the data. While it is obvious that this is indeed 
a gravitational microlensing event, such a poor fit might cast doubt on the extended source 
interpretation. Fits to data of magnification A < 10, which should be negligibly affected by 
extended source effects, display the same excessive scatter around the fit. We therefore conclude 
that the majority of this scatter is contained in baseline measurements, and should not adversely 
affect our interpretation of the data. 

Figure || also presents a scaled schematic of this microlensing event using event parameters 
derived from fits with a limb-darkened source star (fit 3). Included are the lens's Einstein radius 
expressed as t/2, the source radius scaled to 0.0756 Re, and the trajectory of the lens across the 
extended face of the source star. 

Figure |3| provides a visual comparison between the best standard microlensing fit to the data 
(dashed line) and microlensing of a limb-darkened extended source (solid line). Each passband 
was fit independently with a baseline parameter, while event parameters were derived from joint 



fits on the combined dataset. All fits were performed in the MINUIT environment ( James 1994 ). 



5. Determination of Stellar Parameters 

Analysis of the stellar parameters of the source star in MACHO 95-30 (hereinafter referred 
to as the source star) is required if we are to use the additional information provided by u* to 
characterize the lensing object. The most important parameters are its distance (L) and radius 



(i?*). Also valuable are estimates of the source temperature and surface gravity, which will help 
define the limb-darkening coefficients. 

Figure [l] shows the fiducial line for the red giant branch of the metal-rich globular cluster 
NGC 6553 from |Ortolani et al. (199C| ), plotted over the optical CMD for a 2x2 arcmin field 



surrounding the source star. We note that the horizontal branch (HB) in this field is inclined 
like that of NGC 6553, due to the effect of differential reddening. Even though the morphology 
of the red giant branch of this cluster is similar to that of the MACHO field surrounding the 
microlensing event, the magnitude of the source star is fainter than expected, after we account for 
the different reddenings and distances to these sources (we adopt a bulge distance of 8 kpc). The 
smaller apparent luminosity could be due to differences in age, metallicity, or distance. From this 
diagram, we tentatively conclude that this star is either more metal-rich, or more distant than 8 
kpc (or both) than the giants of similar color in NGC 6553. 



5.1. Reddening 

The reddening is patchy in the MACHO field where the source star is located. Figure [| shows 
the amplitude vs. color of all the RR Lyr that we found in this field. The reddening vector is 
horizontal in this figure, so the color spread is mostly due to differential reddening within the field, 



as the RRab sequence in such a diagram should be very tight (Bono et al. f997). There appear 
to be two zones of different obscuration, with E{B — V) ranging from 0.3 to 0.8 mag. The source 
star appears to be located in the region with smaller reddening, as measured from the colors of the 
three closest RR Lyr. We identified a total of two bulge RR Lyr type ab, and one RR Lyr type 
c, within 2 arcmin of the source star. The locations of these 3 RR Lyr are shown in the optical 
CMD (Figure HJ) as triangles. The mean color of the two RRab is V — R = 0.54, and the mean 
color of the RRc is V — R = 0.46. From these colors we deduce E(V — R) = 0.34 ± 0.05. Using the 
extinction law of |Rieke & Lebofski (198^ ), we obtain Ay = 1.36, E(B - V) = 0.44, A K = 0.15, 
E(V -K) = 1.21, E(J -H)= 0.15, E(H - K) = 0.08, and E(J - K) = 0.23. The observed and 
dereddened photometric measurements of the source star are listed in Table |I[ 



5.2. Infrared Photometry 

We observed the field of the source star with the CIRIM camera on the CTIO 1.5m telescope. 
The f/13.5 configuration yields a pixel scale of 0.65 arcsec pixel -1 . The detector is a 256x256 
HgCdTe NICMOS 3 array. Conditions on the night of September 24, 1996 were photometric. 
We used the J, H, and K s filters, and a five-position on-source dithering pattern optimally to 
expose the target field and collect the necessary frames for sky flats. Six 10-second coadds were 
obtained at each position in the dithering pattern for total exposures of 300 seconds in each filter. 
The individual IR images were reduced following standard procedure using dark and sky frames, 



then shifted and combined to create the final images, which cover about 2 arcmin on a side. 
The photometric calibration was obtained by observing the UKIRT Faint Standard Star no. 1 
(G158-100) just prior to observing the field of the source star. The standard star and program 
field were observed at identical airmasses (X = 1.2), thus no extinction corrections were necessary. 
We find K = 9.98, (J — K) = 1.12, and (H — K) = 0.26 mag for the source star. We estimate 
an uncertainty of 0.05 mag in J, H, and K from the photon statistics and aperture corrections. 
Figure |5| shows the optical-IR color-color diagram of stars one magnitude brighter than the 
horizontal branch in the 2x2 arcmin field surrounding the source star, along with the sequences 



corresponding to field giants, field dwarfs, and Baade's Window giants from Frogel &: Whitford 



(1987 ). The colors of the source star are similar to Baade's Window M giants, and we conclude 



that it is unlikely to be a field dwarf. 

Figure || shows an expanded version of the reddening corrected Kq vs (V — K)q optical-IR 
CMD for the bright stars found in this field. Overplotted are the fiducial lines corresponding to 
the red giant branches of the metal-rich globular clusters 47 Tuc from [Frogel et al. (198l| ) and 
the metal-rich globular cluster NGC 6553 from Puarnieri et al. (1996| ), assuming a distance of 
8 kpc for both populations. Figure shows again that the source star is relatively faint with 
respect to the giants in these clusters. This suggests that the star is either more metal-rich 
than [Fe/H] = —0.4, and/or more distant than 8 kpc. As an additional comparison, we have 



plotted data from [Frogel fc Whitford (1987|) for spectral type M3-5 giants from Baade's Window 
(assuming a distance of 8 kpc), and our own estimated fiducial line for these stars. If we assume 
a solar composition for the source star {[Fe/H] = 0), we obtain a distance of about D = 10 kpc. 
The uncertainty in this distance is large; it assumes that all of the scatter about the fiducial line 
of BW giants is due to the line-of-sight depth of the bulge. It is possible that the source star is 
on the metal-rich tail of the distribution of bulge giant metallicities. Comparing the two globular 
cluster tracks in Figure ^ gives one a sense of the dependence of the red giant branch fiducial 
lines on metallicity. If the location of the source star in the CMD were entirely due to metallicity, 
it would be an extremely metal-rich star, which is less probable than a combination of distance 
and metallicity effects at work. Given the position of the source star in the CMD at the outer 
envelope of giants observed in Baade's Window, we conclude that the source star is likely located 
on the "far" side of the bulge. Note that we would expect a giant star in the Sagittarius dwarf 
galaxy to be about 2 magnitudes fainter than the observed Kq magnitude of the source star (see 



also Section 5.3.3J below). 



5.3. Spectroscopy 

We have obtained spectra of the source star in several different observing runs, as will be 
discussed in Section [7[ Here we analyze the spectra taken at the 4m Blanco telescope at CTIO 
with the R-C spectrograph and the Loral 3iCi detector during the nights of 14 and 15 August, 
and 27 September, 1995. The first two spectra were taken 0.5 days before and after, respectively, 



10 



the peak of the microlensing event, while the lens was still in transit across the source face. The 
wavelength coverage of these spectra is from A6230 to A9340 A. A HeNeAr comparison lamp was 
used, which provided the wavelength calibration. The exposures were 60 seconds long, and were 
taken at an airmass of 1.7. All the spectra were reduced in the standard way. The dispersion was 
lA pixel -1 , and the resolution was 5A, as measured from the FWHM of the lines. 

The additional spectrum from 27 September, 1995 was taken at an airmass of 1.1, with a 
total exposure time of 500 seconds. No independent wavelength calibration was obtained for this 
spectrum. The dispersion was 2A pixel -1 , and the resolution was 8A. The wavelength coverage of 
this spectrum, shown in Figure 0, was from A3890 to A9830 A. This spectrum, taken 42.3 days 
after peak magnification, is very similar to the near-peak spectra. This further supports the 
microlensing interpretation, where the overall colors or spectral type must not change. 



5.3.1. Spectral Typing 
We determined the spectral type of the source star by direct comparison with the digital 



spectral atlas of late-type stars of Turnshek et al. (1985 ) and Kirkpatrick et al. (1991 ). This 
comparison was done visually after rebinning our CTIO spectra to lower resolution, and paying 
close attention to the match of the prominent TiO bands. There is a very close match of the 
spectra shown in Figure |?] with a spectral type M4 III, and we adopt this classification for the 
source star. The uncertainty in this classification is estimated to be one spectral subtype. 

Using the spectral type we can check our estimate of the reddening of the source star by 
adopting the unreddened optical and infrared colors of the same type bulge giants. For an M4 
III bulge star, |Frogel fc Whitford (1987J ) give (J - K) = 0.90 and (H - K) = 0.18. and from 
the observed infrared colors we compute E{ J — K) = 0.22, and E{H — K) = 0.08, respectively. 
These reddening values are in excellent agreement with the reddening derived from the optical 
photometry of the RR Lyr in this field (see Table ||) . 



5.3.2. Call Triplet 

The CTIO spectra allows us to determine independently the surface gravity (g) of the source 
star. The spectral region covered includes the Call triplet, which is strong and very sensitive to 
log g. Other spectral lines in the region from A6400 A to A9200 A can be used for this purpose in 
M-type stars, like the Nal doublet at AA8183, 8195A (Kirkpatrick et al. 1991). However, the lines 
of the Call triplet are the most suitable ones, because they have been studied and calibrated by 
a number of authors ( pones et al. 198~4 ; [Humphreys fc Graham 1986 ; Diaz et al. 198£; Alloin & 
Bica 1989| ; |Jorgensen et al. 1992J ; |Erdelyi-Mendez & Barbuy 199 1| ). 



We measured the equivalent width of the Call lines at A8498.06 A, A8542.14 A, and A8662.17 



11 



A in the three spectra obtained with two different resolutions, obtaining EW = 11.0 ±0.5 A, where 
the uncertainty is dominated by the location of the continuum. We used the local continuum 
bands at A8480 A, A8635 A, and A8905 A, defined by pones et al. (1984Q , |Diaz et al. (19891 ), and 



Alloin fc Bica (1989| ), in order to make a direct comparison with their measurements and models. 



Note that the resolution of our spectra is similar to these other works. 

We have also compiled published data (spectral type, T e g , log g, and Call— EW) on M-type 
stars, covering a wide range of temperatures and gravities. In some cases only the two strongest 
lines of the Call triplet are measured, and we have scaled these accordingly in order to include 
the weaker line. Figure |8| shows the equivalent width vs M spectral subtype for dwarfs (V), giants 
(III), and supergiants (I). The position of the source star is shown with the box. This figure 
confirms the conclusions drawn from the IR photometry and spectral classification — the source is 
a giant star. 

Figure [|| shows equivalent width vs log g. The dotted lines show the ±lcr measurements for 
the source star. From this figure we conclude that log g = 1.0 ± 0.2. A similar value is obtained 
by applying the empirical relations found by Alloin fc Bica (1989 ) and Diaz et al. (1989| ), and by 



slightly extrapolating the models of Jorgensen et al. (1992| ) and Erdelyi-Mendez & Barbuy (1991). 



5.3.3. Radial Velocity 

We have also measured the radial velocity of the source star. Unfortunately, we did not 
measure radial velocity standards, and have to rely on a wavelength calibration based on the 
HeNeAr lamps. From the two CTIO spectra of the August run we obtain V r = — 80 ± 5 km s _1 . 
Even though the spectra from the August run have high enough resolution to measure velocities 
good to < 10 km s — 1 , we consider that the most important source of uncertainty is the zero point. 

We consider the possibility that the source star may be located in the Sgr dwarf galaxy. The 



Sgr dwarf galaxy discovered by Ibata et al. (1995 ) is much more extended than previously thought. 



RR Lyr belonging to this galaxy have been identified in fields close to the Galactic plane ( Alard 



et al. 1996| ; Alcock et al. 1997). In fact, the Sgr galaxy extends as far as the field containing the 



source star, as shown by the one Sgr RR Lyrae type ab discovered in this field by Alcock et al. 



(1997). However, the radial velocity measured here rules out membership to the Sgr dwarf galaxy, 



which has V r = +150 km s . 

The radial velocities of bulge giants in off-axis bulge fields have been measured by [Minniti 



(1996 ) and Minniti et al. (1996J) . At the position of the MACHO field, the mean heliocentric 



radial velocity of bulge stars is estimated to be about 25 km s , with a velocity dispersion a = 80 
km s _1 . The radial velocity of the source star is, therefore, consistent with the velocities of bulge 
giants, supporting our assumption that this star belongs to the bulge population. 
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5.4. Stellar Parameters 



We will derive the stellar parameters M, T e g, L, and R for the source star in a variety of 
ways and compare the different values to estimate our uncertainties. When needed, we assume the 
solar values: T & = 5730K, log g® = 4.44, and Mb o i0 = 4.72 mag. The following equations relate 
various stellar parameters: 

(L/L Q ) 

-2 



(R/R Q ) 2 (T eS /T Q )* 



(g/g & ) = (M/M Q )(R/R Q ) 

We begin with the most simple method, adopting the typical stellar parameters of an M4 
giant in the galactic bulge. We would expect the mass of the source star to be ~ IMq, the mass 



appropriate for a bulge giant according to the estimated age of the bulge (e.g. Holtzman et al. 



1993 ). Using only the spectral type information (M4 III), we obtain T e g 
from Lang (1992 ), which gives R = 85Rq. Also, for an M4 III star, R = 



= 3430, M bol = -2.7 
83i? , and T cS = 3600 

from the recent calibrations of Dyck et al. (1996 ). These values are averages over several spectral 
sub-types. The dependence of radius on spectral type is such that for an M3 III star R ss 60R & , 
and for an M5 III star R m 100R e . 



The IR colors allow us to measure the effective temperature. 
T cS = 3900 ± 270 from the calibration of |Feast (19961 ). For V - K 



For J - K Q = 0.89, 
= 5.03, T cS = 3600 from the 
recent calibration of |Bessell et al. (1997|) . The optical photometry also allows us to estimate the 
temperature of the source star. For (V — R)o = 0.93, we find T e g = 3800 ± 200 by differential 
comparison with the giant star 111-17 member of the globular cluster NGC 6553 analyzed by 



Barbuy et al. (1992 ). Comparing these values of temperature with those derived from the 
spectral type information, we adopt a temperature of T e g = 3700 ± 250 for the source star, or 
log(T/T ) = -0.19 ±0.03. 

Numerous bolometric corrections exist in the literature. We find BCk(J — K) = 2.6 and 
BC K {V-K) = 2.7 from |Frogel fc Whitford (1987|) and [Bessell et al. (1997J) , BC H {M giants) = 2.6 
from |Bessell fc Wood (1984]) , and BC V (M4III) = -2.2 ( |Lang 1992| ). These give m hoX = 12.4, 
12.5, 12.6, and 12.6, respectively. We adopt mbol = 12.5 ± 0.1 for the source star. 

In order to compute the absolute magnitudes, we must assume a distance to the star. 
Adopting a distance of 8 kpc, or m — Mq = 14.52, we obtain the absolute magnitudes My = +0.34, 
M K = -4.69, and M bo i = -2.1 mag. This gives log(L/L ) = 2.69 ± 0.04. We can immediately 
derive the radius, \og(R/ 'R®) = 1.72 ± 0.06 or R = 53 ± 8Rq. However, the uncertainty here is 
underestimated given the systematic uncertainty in the adopted distance. From our arguments 



in Section 15^ , we adopt the distance to the source star of 9 ± 1 kpc, giving a distance modulus 
m — Mq = 14.77 ± 0.25. Re-calculating the radius at this distance and incorporating the added 
uncertainty, we find M^ Q \ = —2.25 ± 0.27 and thus R = 61 ± 12Rq. 



Lastly, we derive the mass of the source star using the measured value of log g (Section |5.3| ). 
We find log(g/<70) = —3.44 ± 0.08 and M ~ 0.8 — 2.5M0. The mass derived in this manner is quite 
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uncertain. The largest mass values are unlikely given the age of the bulge t ~ 10 Gyr (Holtzman 



et al. 1993| ). However, a mass slightly larger than 1 M@ is consistent with the source star having 



a high metallicity. 

To summarize, we adopt the following stellar parameters for the source star in MACHO Alert 
95-30, which we categorize as an M4 bulge giant: 

L = 600±200L o , 

T cS = 3700 ± 250K, 

log g = 1.0 ±0.2, 

R = 6l±12R Q , 

D = 9 ± 1 kpc, 

Af « l.OM , 

[Fe/H] « 0. 



6. Determination of Event Parameters 

6.1. Effect of Limb Darkening on the Photometry 

We now integrate an approximate limb-darkening law into the microlensing model, to extract 
more realistic event parameters than those derived with a constant surface brightness disk. 
Limb-darkening coefficients for GMAN's standard R and V passbands exist in the literature for 



T e ff = 3500 K, log g=l for the quadratic form of Equation || ( Claret, Diaz- Cor doves, &: Gimenez 



1995 ; Diaz-Cordoves, Claret, &: Gimenez 1995| ). Limb-darkening coefficients were also calculated 



for MACHO 's non-standard passbands ( plaret 1996 ). These coefficients are listed in Table ||. 



Including this brightness profile in the extended source model further improves the fit x 2 by 9 
(see Table [|, fit 3). While this improvement is formally significant, its interpretation here appears 
unclear. Comparison of fits 2 and 3 in Table || shows no clear trend between data sets, which 
effectively washes out any overall conclusions about the significance of the model. However, we do 
regard the limb-darkened parameters as the more realistic interpretation of the data. Comparison 
of optical and infrared photometry during a lens transit should detect significant color terms as a 
result of source limb-darkening ( Gould fc Welch 1996| ), which would then provide a more robust 
test of stellar atmosphere models. 
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6.2. Implications for Lens Mass 

The limb-darkened fit to finite source microlensing provides a direct measurement of 
u* = 8*/9e = 0.0756, i.e. the angular size of the star as a fraction of the Einstein ring angle. To 
convert this into physical units, we need an estimate of the angular radius 0* of the source star. 
The analysis in Section |5| indicates R* = 61Rq and L = 9kpc, giving 6* = 31.5/iarcsec and thus 
9e = 0.417 milliarcsec. The Einstein diameter crossing time t is measured similarly to the usual 
point-source case, so this provides the proper motion of the lens with respect to the source star 
uj = 26 E /t = 21.5kms _1 kpc -1 = 4.51 mas/yr. 

Thus, an extended source microlensing event provides us with 2 constraints t and to on the 
3 unknowns M, x and v± of the lens. In the following, it is convenient to define v = ujL = v±/x 
to be the relative velocity of the lens projected back to the source plane (cf [Han fc Gould 1996| ); 
substituting te = vxt/2 into Eq. ^, we obtain a one-to-one relationship between the lens mass and 
distance, 

,,. . v^c 2 x , A . 

M ^ = 16GL—x (4) 

(By symmetry, this is the same as Eq. 6 of [Alcock et al. (1995b| ) exchanging observer and source, 
i.e. with v instead of v, and x <-» 1 — x.) For the observed values of v = 193 km s, t = 67.3 
days, we have M{x) = 0.192 M Q x/(l — x), thus the lens may be either a low-mass star roughly 



half-way to the source, or a more massive star closer to the source. Figure 1C shows M(x) for the 
event parameters given in Table || (fit 3), and the radius and distance of the source star obtained 
in Section |5l Since our proper motion error is dominated by uncertainty in the source radius, we 
also include M(x) contours for R* = 61 ± 12R & . 

We may obtain additional constraints on x by using a model for the distributions of sources 
and lenses, since the likelihood of obtaining the observed value of v is also a function of x. 

For given lens mass, the rate of microlensing is proportional to 



dT oc \Jx(l - x) p L {x)v±f s {-v s ) f L {-v L )dxdv s dv L . (5) 

where pi, is the density of lenses at distance x, /l(v&) and fs(ys) are the lens and source 
velocity distribution functions (normalized to unity) in the plane perpendicular to the line of 
sight. The source and lens velocities vg, v^ are related to v by vl = (1 — x)v + x(vs + v), where 
v = (vcos4>,vsiii(j)) and <j) is the (unknown) direction of the relative proper motion. 

Given a model for pL,fs,fL, we may integrate Eq. [B] and thus obtain joint probability 
distributions for any of the variables. Since we have measured v, and the lens mass depends on 
only one unknown (x), we need to consider the joint probability distribution of events in the (x,v) 
plane, and then marginalize to get a probability distribution of x given the observed value of v. 



Thus, we change variables from dv^ to dvd(p, giving 

d(v, 0) 



dT oc \ x(l - x) p L (x) v±f s (v s ) /l(vl) dx dv< 



dvdcfi. (6) 
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We then substitute for v^ and v± and integrate over the unknowns vg and eft, giving a 
likelihood as a function of distance for given v, 



C(x; v) 



ex yJx(l-x)p L (x) v 2 x 3 f f f s (vs) h ((1 - x)v Q + x(v 5 + v)) dv s d(f>, (7) 



dxdv 

This result may be understood as follows: the integrals are over all combinations of source 
and lens velocity which give rise to the observed v. The y/x(l — x) and pl terms arise from 
the x— dependence of the Einstein radius and density of lenses, respectively. There is a factor 
of v± = vx because a given lens contributes a lensing rate ex v±, and a factor of x 2 v from the 
Jacobian det(dvL/d(v, (p)). 

To evaluate Eq. [5], we adopt a disk velocity dispersion of 30kms _1 in each direction, with 
a flat rotation curve of 220 km s -1 . We adopt a bulge velocity dispersion of 80kms -1 in each 
direction, and no bulge rotation. For the density profiles, we use a standard double-exponential 



disk, and a barred bulge as in [Han &; Gould (1995]) , The source is located at galactic coordinates 



I = 3.73°, b = —3.30°, and we assume L = 9kpc, Ro = 8kpc, so the source is behind the galactic 
center. We assume the source is a member of the bulge, so we take fs to be the above velocity 
distribution function of bulge stars, and then we evaluate Eq. |7| separately for p^, fa appropriate 
for lenses either in the disk or the bulge. 

Figure 10 shows the result of Equation [7| as a function of lens distance, for the observed 
v = 193kms _1 . The lower solid line shows disk lenses only, and the upper solid line shows the 
sum of disk + bulge lenses; the relative areas indicate that there is about an 80% probability that 
this lens belongs to the bulge. 

We define a median distance x me d such that half of the integrated likelihood arises from 
x < x me d, giving x me d = 0.77 and a median mass estimate of M me d = 0.67 M©. Similarly, we define 
a 90% confidence interval (x\,X2) such that 5% of the integrated likelihood arises from each of 
x < x\ and x > x% the resulting interval is 0.52 < x < 0.943, which translates into a mass interval 
of 0.21 to 3.2 M©. These results are relatively insensitive to the details of the galactic model, since 
the constraints on x are dominated by the drop in p Bulge f° r x ^ 0.6 and by the geometrical factor 
y/x(l — x) for x — > 1. 



For reference, Figure 11 shows v J C(x;v) dx, i.e. the lensing rate per unit logu for disk and 
bulge lenses. The observed value of v = 193kms _1 is well within the range of expected values. 
The areas under the two curves reflect the fact that bulge lenses produce about 3 times the event 
rate of disk lenses for the above models. This figure shows that there is a large overlap in the 
distributions; events with v <^ 150kms _1 are produced almost entirely by bulge lenses, but the 
more common events with larger v arise from both disk and bulge lenses; thus, proper motion 
measurements towards the bulge are not as useful as parallax measurements for constraining the 



lens distance. (This has been previously noted by Han &: Gould (1995|) ) 



Note that this is distinctly different from the LMC case, where either proper motion or 
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parallax measurements provide a good separation of the various lensing populations (Gould 1996) 



6.3. Constraints on Lens Luminosity 

As seen above, if the lens is more distant than x > 0.9, it has a relatively large mass. If it is a 
hydrogen-burning star, we can constrain this possibility as follows. For the above proper motion 
of uj ~ 4.5 mas/year, the lens is completely unresolved from the source during our observations, 
thus the observed light-curve of our 'object' is simply the sum of that from the source, the lens 
and possibly other superposed stars. 

Any flux from the lens would add a constant un-magnified baseline which would distort the 
microlensing fit. Although in most cases this additional flux would be a small fraction of the flux 
of the source star, with the high precision measurements here this would be detectable even if the 
lens is considerably fainter than the source. 

We have fit "blended" microlensing to data of magnification A < 10. This fit is of the form 
fi(t) = A(t)foi + fui, for i = each passband, where /o« is the baseline flux of the source star and 
fui is the total flux of unlensed stars superposed on the source, i.e. the sum of flux from the lens 
and any other superposed stars. 

We note that the Xdof °f an light-curves in these blended fits are similar to those in the 
extended source fits. We therefore exclude the MACHO data in our lens brightness estimate due 
to excessive baseline scatter; however, this turns out to be unimportant. The most follow-up data 
with A < 10 are contained in the CTIO and UTSO observations. For the CTIO R, UTSO R, and 
UTSO V passbands, we find fm/foi « (-0.021 ± 0.016), (-0.019 ± 0.017), and (-0.015 ± 0.016), 
respectively, consistent with zero flux from the lens. For completeness, we note that similar results 
are obtained from the MACHO data, with fm/foi ~ (-0.012 ± 0.017) and (-0.026 ± 0.018) for 
MACHO V and MACHO R, respectively. MJUO and Wise data did not contain enough baseline 
observations to determine accurately the amount of unlensed flux. In the following, we take 
± 0.02 as a conservative limit. 

If the lens is a main-sequence star, we may predict its apparent brightness using the 
mass-distance relation of Eq. |4J. Assuming a main-sequence V-band mass-luminosity relation 
L/L & = (M/M ) 3 ' 5 , and 1.5 magnitudes of extinction, we find that the apparent brightness 
of a main-sequence lens would be e.g. 1% of that of the source star for x = 0.83, and 10% for 
x = 0.90. Note that the implied lens brightness increases with distance since the rapid rise in 
Mix) outweighs the r~ 2 term; note also that a giant lens is excluded at any distance. 

Instead of applying a sharp limit x < 0.88, it is more rigorous to multiply the likelihood 
function of Eq. by the Gaussian probability that the lens brightness is consistent with the above 
constraint on fu/fo- This causes a rapid roll-off in the likelihood function for Dl > 7.5 kpc, as 
shown by the dotted curves in Figure [jij. Assuming a main-sequence lens, we can re-compute the 
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constraints on x and M from the likelihood function with the brightness factor, giving a median 
x = 0.73, a mass of M me( i = 0.53 M , and 90% confidence intervals of 0.48 < x < 0.84 and 
0.18 < M < 1.05 Mq. Thus, the inclusion of the lens brightness constraint reduces the median 
mass only slightly, but considerably strengthens the upper limit. 



7. Spectral Variation during MACHO 95-30 

In addition to the CTIO observations discussed in Section [T^, we have obtained spectra at 
Mount Stromlo and Keck [] Observatories. Table || lists the complete catalogue of spectroscopic 



observations of this event. Figure 12 schematically shows the location of the lens with respect 
to the source when these different observations were made. Nightly spectra of the source star 
were taken between August 18 and 25, 1995 (from 2.8 to 9.8 days after peak magnification), 
with the Cassegrain spectrograph at the MSO 74-inch telescope. They cover the wavelength 
interval 6240-6770 A, with a dispersion of 0.90 A per pixel, and a resolution FWHM = 4.6 A as 
determined from the comparison lamp spectra. The total exposures were typically 1000 seconds 
long. These MSO spectra are wavelength calibrated and sky subtracted, but not fluxed. The S/N 
per resolution element is listed in Table ^, along with other relevant data. 

The CTIO and MSO spectra have similar resolution, with the exception of the September 27 
CTIO spectrum, which has lower resolution. We rebinned the CTIO and MSO spectra in order 



to compare them directly. Figure IS shows the spectral sequence, which constitutes one of the 
most extensive and homogeneous series of spectroscopic observations of a microlensing event to 
date (see also Benetti et al. 1995|) . It is clear that the microlensing does not change the spectral 



type of the star, nor does it strongly affect the major spectral features. However, subtle effects 
may appear in some spectral lines, as discussed by Loeb fc Sasselov (1995| ), which warrants a more 
careful comparison among the spectra. 

Before this discussion, we note that the MSO spectrum of August 18 (the first spectra taken 
after the lens transit) is anomalous, showing a dip at A6520A that looks like an unidentified 
bandhead. Although we have examined several possible sources for this feature, including checking 
the flat fields and comparing the spectra of other stars in the field, the cause of this dip is still 
unexplained. 

Each spectrum was divided by the median combination of all spectra (this operation was 
repeated before and after continuum subtraction, in order to account for possible differences due 



to flux calibration of the CTIO spectra but not of the MSO ones). Figure 14 shows the sequence 
of these ratios. While there are no strong deviations, the two spectra taken closest to maximum 
magnification show stronger Ha than the rest. Also, these two spectra show stronger TiO bands 
than the following ones. The total TiO absorption, however, started to climb steadily after 18 



5 The Keck telescope project was made possible by a generous grant from the W.M. Keck Foundation. 



18 



August, 1995. We measured the equivalent widths of several spectral features of interest within 
6240-6770 A with the SPLOT package within IRAF. Figure [0] shows the equivalent width of Ha, 
and the combined intensity of 4 TiO bands in the 6650-6750 A region. Included are conservative 
error bars based upon uncertainty in placing the local continuum. The equivalent width of the Ha 
line was scaled to the equivalent width measured in the HIRES spectra of 18 and 19 August, 1995. 
These equivalent widths were measured with respect to the local pseudo continuum, in the same 
way for all the spectra, in order to avoid systematic effects. 

While we cannot explain in detail the behavior of these spectral features, changes in the 
equivalent width of these lines have been predicted by Loeb fc Sasselov (1995| ). These changes are 



due to limb brightening effects in the cores of resonance lines (like Ha) due to a very extended 
photosphere. Resonant line scattering may also affect the optical TiO bands in the case of M 
giants like the MACHO 95-30 source ( [Basselov 1997| ). Detailed modeling of the spectra of this 
event in particular is needed, as this will help with the interpretation of forthcoming similar events. 

Three spectra were also taken in the nights of August 18, August 19, and August 21, about 
5 hours later than the MSO spectra, using HIRES at the Keck 10-m telescope. These high 
resolution spectra also have high S/N, as listed in Table pi which includes the S/N per resolution 
element of the lowest and highest orders. The spectra of August 18, August 19, and August 21 
consist of 30, 25, and 40 orders, respectively. The blue region below about 4000A covered in the 
August 21 spectrum does not have much information due to low counts. We have checked different 
resonance lines (Call, Ha, H/3, H7, etc), finding no large variations. The August 21 spectrum 
has coverage extending to the blue, including the Ca H and K lines, but misses the Ha line. The 
spectra of August 18 and 19 do not include these Ca lines. The Ha equivalent width differs by 
only 4% between August 18 and 19, as measured from the high quality HIRES spectra (1.087 vs 
1.039A). This variation occurs in the core of the line. Otherwise, there is no significant change 
between the Ha profiles of August 18 and 19. The strong TiO band heads at AA 4954, 5166, and 
5447 A are present in all three HIRES spectra. The difference of about 10% in strength of these 
bands between August 18 and 21 seen in the MSO spectra is confirmed. 

We have also fit the Nal (D) line profiles in the HIRES spectra with an interstellar cloud 
model which uses 12 clouds over the velocity range — 60. 8 km s -1 to +202.3 km s . The clouds 
with extreme velocities only fit one line of the Na doublet because of problems with the continuum 
determination and blending with other lines. The interstellar clouds indicate that the star is at a 



great distance and therefore it is certainly a giant, as discussed in Section 5.3 



We also measured the radial velocities of the MSO spectra using the cross-correlation routine 
FXCOR within IRAF. From the MSO spectra, we measure a mean velocity of V r = — 76±4kms~ 1 , 
which is in excellent agreement with the mean velocity of V r = — 80 ± 5kms _1 from the CTIO 
spectra, reduced and calibrated independently (Section ^). The velocities of the eight MSO 
spectra agree with each other within an rms of 12kms _1 . We measured also the radial velocity of 
the August 21 HIRES spectrum, V r = —86.2 ± 3.9kms _1 , using the strong Cr lines at AA 4254, 
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4274, and 4289 A. All these velocities are in excellent agreement, and we adopt a final value of 
V r = -81 ± 5kms _1 for the source star in MACHO 95-30. 

Subtle radial velocity variations (of the order of 1 — 2kms _1 ) may also be expected during 
microlensing events such as this one ( Maoz fc Gould 1994J ), and can give an independent 



measurement of the proper motion. Alternatively, combining the radial velocity variation with the 
extended source magnification effect, one can determine the projected rotation velocity of the star 



(|Gould 1997|) . We note that detecting such velocity differences may be possible using the HIRES 



spectra, after detailed corrections and modeling that are beyond the scope of the present paper. 

In summary, we have obtained a large number of spectra of the microlensed source star during 
MACHO Alert 95-30. These spectra show subtle variations, which support the predictions of 
existing microlensing and stellar models. These spectra demonstrate that differential magnification 
of an extended source star during a gravitational microlensing event may be used to probe 
otherwise inaccessible stellar atmosphere "fine structure". 



8. Conclusions 

We have observed the breakdown of the point-source approximation for the first time in a 
"normal", single lens microlensing event. This deviation was anticipated based upon real-time 
information provided by MACHO-GMAN, and dynamically scheduled for GMAN follow-up 
observations throughout the event. In effect, we have "resolved" a star which is ~ 9kpc from 
Earth. The resulting deviation from point-source microlensing provides a second constraint 
equation for the system, partially reducing the ambiguity between lens mass, position, and 
transverse velocity. 

Measurement of the effective radius of the source star during the microlensing event provides 
an estimate of the scale of the lens's Einstein radius. With additional information about the radius 
and distance of the source, we determine the angular velocity of the lens with respect to our line 
of sight to the source. This allows a statistical estimate of the lens mass by assuming velocity 
distributions of sources and lenses. We conclude that this event is due to lensing by an object of 
mass O.67^q' 46 M0 (at the 90 % confidence level), such as a white dwarf or a neutron star, which 
is most likely in the bulge (80 % c.l.). If the lens is a main-sequence star, the upper end of this 
mass range is excluded, and our mass estimate becomes O.53^ O ' 35 M0. The results of this paper 
are attributable to the coordinated efforts of GMAN observers, and are a robust endorsement for 
microlensing follow-up implemented at 1-m class telescopes. 

We have also detected for the first time subtle variations in the spectra of a star undergoing 
microlensing. Anomalously strong equivalent widths of Ha and TiO were detected while the lens 
was in transit across the source, an effect anticipated due to varying line widths across the face of 
the star. Given our ability to predict limb crossings in real-time, this technique may be used to 
map out the center-to-limb variations in the spectrum of a future microlensed source. 
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Fig. 1. — Single epoch optical V vs V — R color-magnitude diagram of the field surrounding the 
source star, which is indicated with a filled circle. Three RR Lyr used to estimate the reddening 
are plotted as triangles. The fiducial loci of the giant branch of the metal-rich globular cluster 
NGC 6553 {[Fe/H] = -0.2) from Ortolani et al. (1990) is also indicated. 

Fig. 2.— 1995 light-curve of MACHO Object 101.21821.128 (Alert 95-30). The data represent 
MACHO-GMAN observations of this event in all passbands within a window of 70 days around 
peak magnification. The date this event was detected by the MACHO Alert system is indicated 
with an arrow. An additional schematic relates the scale of the lens's Einstein radius to the angular 
size of the source star, and indicates transit of the lens across the source face. 

Fig. 3.—- Peak of MACHO Alert 95-30, showing the best standard microlensing fit to the data 
( ) and an extended source microlensing fit incorporating source limb-darkening (— — ). 

Fig. 4. — Amplitude vs. color of the RR Lyr type ab in the MACHO field where the event is 
located. The two RRab closest to the source star are indicated with asterisks. We find Ay = 1.35 
mag. 

Fig. 5. — Optical-infrared (V — K )o vs (J — K)q color-color diagram for the giants brighter than 1 
mag above the horizontal branch in the field surrounding the source star. The fiducial loci of field 
giants, dwarfs, and Baade's Window giants from Frogel & Whitford (1987) are indicated with the 
solid, short-dash, and long-dash lines, respectively. 

Fig. 6. — Optical-IR Kq vs (V — K)q color-magnitude diagram of the field surrounding the 
source star. The fiducial loci of the giant branches of the metal-rich globular clusters 47 Tuc 
{[Fe/H] = -0.7), and NGC 6553 ([Fe/H] = -0.2) from Guarnieri et al. (1996) are indicated. 
Additionally, we have plotted BW M3-5 giants (Frogel & Whitford, 1987) as small crosses and 
drawn our own estimated fiducial line through these points. 

Fig. 7. — Baseline spectrum of the source star obtained at CTIO on September 27, 1995, 42 days 
after maximum magnification. The flux calibration is unreliable beyond 9000A. 

Fig. 8. — Call equivalent width in A vs M spectral sub-type for dwarfs (circles), giants (asterisks), 
and supergiants (triangles). The box encloses the drier measurements of the source star in MACHO 
95-30. 

Fig. 9. — Call equivalent width in A vs log g, with the dotted lines indicating the Call equivalent 
width of the source star in MACHO 95-30 (drier). 

Fig. 10. — Lens mass plotted as a function of its distance, using Equation EL This assumes a 61i?0 
source, with error contours provided at ±12Rq, and event parameters listed in Table 13] (fit 3). The 
solid likelihood curves show relative probabilities for disk and bulge lenses as a function of their 
distance. The dotted lines include an upper limit on the brightness of a main-sequence lens. 
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Fig. 11. — Lensing rates per unit log v for disk and bulge lenses, indicated with the solid and 
dashed lines, respectively. 

Fig. 12. — Spectral observations listed in Table y plotted on the source plane. The relative sizes 
of the star (solid circle) and the lens's Einstein radius (dotted circle) are plotted to scale, in units of 
Re- The solid points show the position of the lens at the different times when spectra were taken. 
The observatories where the spectra were taken are indicated with the tick marks. 

Fig. 13. — Spectra of the source star taken at CTIO and MSO. All of these have similar resolution, 
except for the top one. Note Ha at A6553 A, and the TiO band-heads at AA 6647, 6676, 6711, and 
6742 A. Days from the peak of the microlensing event are indicated along the vertical axis. 

Fig. 14. — Same as Figure |13|, showing the spectra divided by the median combination of all. 
These have been displaced in the vertical direction and arbitrarily scaled, such that the flattest 
spectra here actually display the largest range of excess from the median combination. Note the 
changes in Ha and the TiO bands in the two spectra nearest peak magnification. 

Fig. 15. — Temporal evolution of the equivalent width in A of Ha and the TiO band- heads, 
measured with respect to the local pseudo-continuum. The EW plotted in the right panel are the 
sum of the TiO bands at AA 6647, 6676, 6711, and 6742 A. 
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Table 1. Global Microlensing Alert Network 



Observatory 



Latitude, Longitude 



Aperture 



Detector 



Plate Scale ("/pixel) 



MSO a 149 ? 00.5 ™ E , 35 ? 19.2 ™ S 1.27 m 

CTIO b 70 ? 48.9 ™ W , 30 ? 30.9 ™ S 0.9 m 

MJUO c 170 ? 27.9 ™ E , 43 ? 59.2 ™ S 0.61 m 

UTSO d 70 ? 42.0 ™ W , 29 ? 00.5 ™ S 0.61 m 

WISE e 34 ? 45.8 ™ E , 30 ? 35.8 ™ N 1.0 m 



8 x 2048 x 2048 Loral 
1024 x 1024 Tcktronics 
1536 x 1024 KAF1600 
512 x 512 Photometries 
1024 x 1024 Tektronics 



0.63 
0.40 
0.23 
0.45 
0.70 



a Mount Stromlo Observatory, Canberra, Australia. 

b Cerro Tololo Inter-American Observatory, Cerro Tololo, Chile. 

c Mount John University Observatory, Lake Tekapo, New Zealand. 

d University of Toronto Southern Observatory, Las Campanas, Chile. 

c Wise Observatory, Mitzpe Ramon, Israel. 

Summary of observatories participating in CM AN follow— up observations. 
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Table 2. MACHO Alert 95-30 microlensing statistics 



Fit a t b t C U m in U* UJ d 



1 


1321.2 (1) 


67.70 (13) 


0.04069 (12) 





??? 


2 


1321.2 (1) 


67.36 (1) 


0.05579 (1) 


0.07335 (1) 


22.1 (50) 


3 


1321.2 (1) 


67.28 (27) 


0.05408 (20) 


0.07561 (9) 


21.5 (49) 



a Fit 1 is the best standard microlensing fit to the data, Fit 2 
incorporates the extended size of the source star, and Fit 3 takes into 
account limb-darkening of the source. 

b JD - 2448623.50. 

c Einstein diameter crossing time. 

Lens angular velocity, relative to the source, is in kms^ 1 kpc -1 and 
assumes a 61 ± 12Rq source at 9 kpc. 

Comparison of event parameters between microlensing fits for 
MACHO Alert 95-30. Statistics are derived from simultaneous fits on 
all passbands. Reported uncertainties in the final significant digit(s) are 
the maximum extent of the surface in parameter space which has a x 2 
greater than the best-fit value by 1. 
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Table 3. MACHO Alert 95-30 microlensing statistics 



Passband # Observations dm a Fit 1 \ 2 Fit 2 x 2 Fit 3 x 2 



dACHO R 


205 


0.018 


980.33 


847.49 


839.15 


dACHO V 


292 


0.023 


1069.65 


964.16 


963.86 


CTIO R 


103 


0.019 


307.58 


124.97 


126.51 


MJUO R 


41 


0.019 


184.56 


46.68 


42.62 


UTSO V 


55 


0.020 


101.34 


48.64 


50.28 


UTSOR 


96 


0.013 


488.20 


82.91 


90.09 


WISER 


12 


0.013 


80.71 


13.76 


7.22 



TOTAL 804 3212.40 2128.64 2119.76 

a Average error, in magnitudes, for each passband. 

Individual microlensing statistics for MACHO and GMAN observations 
of Alert 95-30. The number of constraints per passband for fits 1,2, and 3 
are four, five, and five, respectively. 
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Table 4. Photometry of the source star in MACHO 95-30 



Observed 


Extinction 


Dcrcddcncd 


Abs Mag, 8 kpc 


Abs Mag, 9 kpc 


V = 16.21 


Ay = 1.35 


Vo = 14.86 


M v = +0.34 


M v = +0.59 


K = 9.98 


A K =0.15 


ATo = 9.83 


M K = -4.69 


M K = -4.45 


V -R= 1.39 


E(V -R) = 0.34 


V - # = 1.05 






J - K = 1.12 


E(J-K) = 0.23 


J - X = 0.89 






H - K = 0.26 


E(H - K) = 0.08 


H - K = 0.18 






V - K = 6.23 


E(V - K) = 1.21 


V - ATo = 5.03 






Bolometric 




BC*k = -2.7 ±0.1 


M bol = -2.0 


M bol = -2.25 
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Table 5. Limb-darkening coefficients for the source star 



Passband 



MACHO V 1.140 -0.284 

MACHO R 0.825 -0.051 

Standard V 1.072 -0.228 

Standard R 0.910 -0.108 

Limb darkening coefficients 
for Equation |3] used 
to approximate the brightness 
profile of the source star in 
MACHO Alert 95-30. 
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Table 6. Summary of spectral observations appearing 


in Figure 12 




JD 


Telescope & Instrument 


Dispersion 


Resolution 


Coverage 


S/N 


2449944.234 


CTIO 4m/RC Sp. 


1.0 A pix" 1 


4.0A 


6230-9340A 


-100 


2449945.211 


CTIO 4m/RC Sp. 


1.0 A pix" 1 


4.0A 


6230-9340A 


-100 


2449947.221 


Keck lOm/HIRES 


0.04 A pix" 1 


0.2A 


4309-6739A 


92-256 


2449947.524 


MSO 74" /Cass. Sp. 


0.9 A pix- 1 


4.6A 


6240-6770A 


106 


2449948.223 


Keck lOm/HIRES 


0.04 A pix" 1 


0.2A 


4835-7282A 


166-296 


2449948.504 


MSO 74" /Cass. Sp. 


0.9 A pix" 1 


4.6A 


6240-6770A 


117 


2449949.509 


MSO 74" /Cass. Sp. 


0.9 A pix" 1 


4.6A 


6240-6770A 


113 


2449950.322 


Keck lOm/HIRES 


0.04 A pix" 1 


0.2A 


3750-6065A 


20-141 


2449950.523 


MSO 74" /Cass. Sp. 


0.9 A pix" 1 


4.6A 


6240-6770A 


91 


2449951.498 


MSO 74" /Cass. Sp. 


0.9 A pix" 1 


4.6A 


6240-6770A 


96 


2449952.508 


MSO 74" /Cass. Sp. 


0.9 A pix" 1 


4.6A 


6240-6770A 


72 


2449953.494 


MSO 74" /Cass. Sp. 


0.9 A pix" 1 


4.6A 


6240-6770A 


78 


2449954.513 


MSO 74" /Cass. Sp. 


0.9 A pix" 1 


4.6A 


6240-6770A 


66 


2449987.000 


CTIO 4m/RC Sp. 


2.0 A pix" 1 


8.0A 


3890-9830A 
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